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We discuss the results of calculations of the shear viscosity and thermal conductivity of pure 
neutron matter, carried out within the Landau- Abrikosov-Khalatnikov formalism. The probability 
of neutron-neutron collisions in the nuclear medium has been obtained from a realistic potential, 
using both the correlated basis function and the G-matrix approach. The results of our work 
indicate that medium modifications of nucleon-nucleon scattering are large, their inclusion leading 
to a dramatic enhancement of the transport coefficients. On the other hand, the results obtained 
from the two theoretical schemes appear to be in fairly good agreement. 



PACS numbers: 24.10.Cn,25.30.Fj,61.12.Bt 



I. INTRODUCTION 

The knowledge of transport properties of neutron mat- 
ter is relevant to the understanding of a variety of neutron 
star properties. Viscosity plays a crucial role in determin- 
ing the onset of the gravitational- wave driven instability, 
associated with the excitation of r-modes, in rapidly ro- 
tating stars [l[ , while thermal conductivity is one of the 
factors determining neutron star cooling [2|. 

Unlike the equation of state (EOS), which is generally 
obtained from realistic dynamical models, strongly con- 
strained from nuclear systematics and nucleon-nucleon 
scattering data, the non-equilibrium properties of neu- 
tron star matter are often studied using oversimplified 
models of the nucleon-nucleon (NN) interaction. 

The main difficulty involved in the calculation of the 
transport coefficients within the formalism originally de- 
veloped by Abrikosov and Khalatnikov Q , based on Lan- 
dau theory of normal Fermi liquids 0], is the determi- 
nation of the nucleon-nucleon (NN) collision probabil- 
ity in the nuclear medium. Most studies of the trans- 
port properties of neutron star matter have circumvented 
this problem neglecting medium modifications of the NN 
cross sections altogether, and using the measured NN 
scattering phase shifts to construct the collision proba- 
bility (EE 0i- 

Nuclear many body theory provides a consistent frame- 
work to obtain the in-medium NN cross section and the 
transport coefficients of nuclear matter from realistic NN 
potentials, using either the G-matrix [8[ or the correlated 
basis function (CBF) [3 formalism. In both approaches 
one can define a well behaved effective interaction, suit- 
able for use in standard perturbation theory in the Fermi 
gas basis and allowing for a consistent treatment of equi- 
librium and non-equilibrium properties [1, 0, [l(| . 

In this paper we discuss the results of calculations of 



the shear viscosity and thermal conductivity of pure neu- 
tron matter, carried out using the CBF and G-matrix 
effective interactions. 

In Section [ill after outlining the elements of nuclear 
many body theory, we analyze the main features of the 
CBF and G-matrix effective interactions, while Section 
Mil is devoted to the discussion of the in-medium NN 
cross section in the kincmatical setup relevant to the cal- 
culation of the transport coefficients. The main features 
of the Abrikosov-Khalatnikov formalism are reviewed in 
Section IIV1 where we also present the results of numeri- 
cal calculations. Finally, in Section [V] we summarize our 
findings and state the conclusions. 



II. EFFECTIVE INTERACTIONS IN NUCLEAR 
MANY BODY THEORY 



Nuclear many body theory (NMBT) is based on the 
tenet that nuclei can be described in terms of point like 
nucleons, whose dynamics are dictated by the hamilto- 



H 



E 



2m 



j>i 



k>j>i 



(1) 



kj and m being the momentum of the i-th nucleon and 
its mass, respectively. 

The nucleon-nucleon (NN) potential Vij reduces to the 
Yukawa one-pion exchange potential at large distances, 
while its behavior at short and intermediate range is de- 
termined by a fit of deuteron properties and NN scatter- 
ing phase shifts. The state-of-the-art NN parametriza- 
tion referred to as Argonne v%s potential [llfl is written 
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in the form 

18 

% = 5>™( r «)°s • (2) 

n=l 

In the above equation 

0^ 8 = [l,(«7 < -«7 i ),Si J -]®[l,(T < -T j )] (3) 

where cr^ and Tj are Pauli matrices acting in spin and 
isospin space, respectively, and 

3 

S H = — (°"' ' V v)( (T 3 " r y") - (°» ' °j) • ( 4 ) 

The operators corresponding to n = 7, . . . , 14 are associ- 
ated with the non-static components of the NN interac- 
tion, while those corresponding to n = 15, . . . , 18 account 
for small charge symmetry violations. Being fit to the full 
Nijmcgen phase shifts data base, as well as to low energy 
scattering parameters and deuteron properties, the Ar- 
gonne Via potential provides an accurate description of 
the scattering data by construction. 

The three-nucleon potential Vijk, whose inclusion is 
needed to reproduce the observed binding energies of the 
three-nucleon system and the empirical nuclear matter 
equilibrium properties, consists of the Fujita-Miyazawa 
two-pion exchange potential supplemented by a purely 
phenomenological repulsive contribution [l2T |. 

The predictive power of the dynamical model based on 
the hamiltonian of Eq. (JTJ) has been extensively tested by 
computing the energies of the ground and low-lying ex- 
cited states of nuclei with A < 12. The results of these 
studies, in which the many body Schrodinger equation 
is solved exactly using stochastic methods, turn out to 
be in excellent agreement with experimental data [HI]. 
Accurate calculations can also be carried out for uni- 
form nuclear matter, exploiting translational invariancc 
and using the stochastic method [l4| , the variational ap- 
proach [l5| , or G-matrix perturbation theory [l6|, [l7| ■ 

One of the most prominent features of the NN potential 
is the strongly repulsive core, whose cleanest manifesta- 
tion is the observed saturation of nuclear charge densities. 
Due to the presence of the core, the NN potential cannot 
be used to carry out ab initio microscopic calculations 
of nuclear observables using standard perturbation the- 
ory. The matrix elements of the interaction hamiltonian 
between eigenstatcs of the noninteracting system, Fermi 
gas states in the case of uniform nuclear matter, turn out 
to be very large, or even divergent. 

In the G-matrix approach the above problem is circum- 
vented replacing the bare NN potential with the well be- 
haved operator G, defined through the Bethc-Goldstonc 
equation 

(ij\G(E)\kl) =G ijM {E) 

= Vij ki + y^v l j, mn — — ^ — — —G mn ki(E),(S) 



where i = (ki,Sj,£i), k^, Si and ti being the momentum 
and the spin and isospin projections specifying the i-th 
single particle state. The Pauli operator Q mn restrict the 
sum over intermediate states to those compatible with 
the exclusion principle, while the so-called starting en- 
ergy E corresponds to the sum of the non-relativistic 
energies of the interacting nucleons. 

The single-particle energy of a nucleon in the state i is 
given by 

* = + Re[Ui] , (6) 

where Ui describes the mean field felt by the nucleon due 
to its interactions with the other particles of the medium. 
In the so-called Brueckner-Hartree-Fock approximation, 
Ui is calculated in the "on-shell approximation" through 
a self-consistent process. The resulting expression is 

Ui= Yl (ij\G(E = e i + e j )\ij) a , (7) 

where the sum runs over all occupied states in the 
Fermi sea {F} and the two-nucleon matrix elements 
are properly antisymmctrized. We note here that the 
so-called continuous prescription [I?) has been adopted 
for the single-particle potential when solving the Bcthe- 
Goldstone equation. As shown in Ref.[l8[, the contribu- 
tion to the energy per particle from three-hole line dia- 
grams is minimized by this prescription. 

Once a self-consistent solution of the G-matrix is 
achieved, the energy per particle at the two-hole line level 
takes the form 

J = l^ + l E (^ = e i + e j )\i j)a , (8) 

where hp is the Fermi momentum, related to the density 
through the relation p = vkp / 67r 2 , v being the spin- 
isospin degeneracy of the momentum eigenstates (y = 
2 and 4 for pure neutron matter and symmetric nuclear 
matter, respectively). 

In the approach based on correlated wave functions one 
uses the bare potential, v, whose non perturbative effects 
are incorporated in the basis states, obtained from the 
Fermi gas states \npo) through the transformation 

|n) = i> FG ) . (9) 

The operator F, embodying the correlation structure in- 
duced by the NN interaction, is written in the form 

F = Sl[f ij , (10) 

ij 

where S is the symmetrization operator accounting for 
the fact that, in general, [fij,fik] 0- The two-body 
correlation functions fy, whose operator structure re- 
flects the complexity of the NN potential, is written in 
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the form 



k^ and k 2 . The total energy of the initial state 
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with the Ofj given by Eq.©. 

Within the correlated basis functions (CBF) approach, 
at two-body cluster level one finds [§, 03] 



(12) 
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and the derivatives act on the relative coordinates. 

One would be tempted to exploit the analogies between 
Eqs.© and (fT2"| . to establish a direct link between G and 
Vcff- Hovcwer, determining such a connection at operator 
level is not trivial. To see this, just consider that, while 
all matrix elements of G involve the bare interaction v 
only, the matrix elements of V e g also include purely ki- 
netic contributions, not containing v, which arise from 
the derivatives of the correlation functions. In addition, 
unlike V e g, G exhibits an explicit energy dependence. 

The CBF effective interaction, being defined through 
its ground state expectation value, is somewhat limited 
in scope, with respect to the G-matrix effective inter- 
action. However, a systematic comparison between the 
two formalisms can be carried out at the level of ma- 
trix elements. In this work we will focus on the ma- 
trix elements of the effective interactions in momentum 
space relevant to the calculation of the NN scattering 
rate in pure neutron matter, whose knowledge is required 
to obtain the transport coefficients within the Landau- 
Abrikosov-Khalatnikov formalism. 

We have used the truncated version of the Argonne 
vis potential referred to as v' 6 [l|| , whose definition only 
involves the static contributions, i.e. those correspond- 
ing to n < 6, in Eq.Q. The CBF effective interaction 
derived from this potential has been also used to obtain 
weak response of nuclear matter at moderate momentum 
transfer [§, [2(| ■ 

For the sake of simplicity, in this work we have ne- 
glected the contribution of the three nucleon potential 
appearing in Eq.([l]). 



III. NN SCATTERING IN THE NUCLEAR 
MEDIUM 

A. Kinematics 

Consider the process in which two nucleons carrying 
momenta ki and k2 scatter to final states of momenta 
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can be conveniently rewritten in terms of the center 
of mass and relative momenta, K = ki + k2 and 
k = (ki k 2 )/2, as 
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(15) 



with M = 2m and [i = m/2. 

In the reference frame in which the center of mass of 
the system is at rest (CM frame) E = Ecm = £re\, while 
in the lab (L) frame, in which k 2 = 0, E = E^ = 2£ Ie \. 

The analysis of the NN scattering rates relevant to the 
calculation of the transport coefficients is carried out in 
the frame in which the Fermi sphere is at rest, often re- 
ferred to as Abrikosov-Khalatnikov (AK) frame. More- 
over, in the low-temperature regime, in which the results 
of Ref. [|[ are applicable, scattering processes can only 
involve nucleons with momenta close to the Fermi mo- 
mentum. Therefore, one can set 

1^1 = 1^1 = 1^1 = 1^1=^. (16) 

At energies below pion production threshold the scatter- 
ing process is elastic, so that the requirement of energy 
conservation 
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2fc|(l + cos0) 

(k' 1 +k 2 ) 2 = 2A;|(l + cos0 / ) , (17) 



implies that the angle between the momenta of the two 
nucleons is the same before and after the collision. In 
general, however, the angle <j> between the initial and 
final relative momenta, k and k' = (k^ — k 2 )/2, defined 
through 
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(18) 



docs not vanish. Hence, for any given Fermi momentum, 
i.e. for any given matter density, the scattering process 
in the AK frame is specified by the center of mass energy 



2m 



(1 



(19) 



and the two angles 9 and 4>. 

As the NN scattering cross section is often evaluated in 
the CM frame, it is convenient to establish a relationship 
between kinematical variables in the CM and AK frames. 
Exploiting the frame invariance of the relative energy we 
easily obtain 
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while the center of mass scattering angle 6cm can be 
indentified with cf>, defined in Eq. (Ti"5|) . 
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B. Cross section 

In both the G-matrix and CBF efffective interaction 
approaches, the NN cross section in matter at density p, 
can be written in the form 
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where m* is the nucleon effective mass, and the transition 
amplitude in the channel of total spin S and initial and 
final spin projections M and M', A4g iM (9, 4"), involves 
the matrix elements of either G or V e g between Fermi 
gas states. 

Numerical calculations of the cross sections are carried 
out expanding AAg [M (9, </>) in partial waves. In the case 
of pure neutron matter, in which the total isospin of the 
interacting pair is T = 1, the expansion only involves 
partial waves of even (odd) angular momentum in spin 
singlet (triplet) states. 

The matrix elements of the CBF effective interaction 
can be written in the form 



M//(fc) 



r 2 drMkr){£'SJ\V eS \lSJ}j e ,(kr) . (22) 



In the above equation, £ and J denote orbital and to- 
tal angular momentum, respectively, ji is the spherical 
Bessel function, \iSJ) is the spin-angle state and r is the 
magnitude of the relative distance. Note that, due to the 
presence of the tensor operator of Eq. ((4|) , the NN poten- 
tial couples states of different orbital angular momentum. 
The above matrix elements are directly comparable with 
those obtained from the partial wave expansion of the 
G-matrix (see Eq.©). 

In Fig. Q]we show the matrix element M( 1 5o) = Mq§, 
evaluated at nuclear matter equilibrium density, po = 
0.16 fm -3 (top panel), and 2po (bottom panel), in the 
kinematical setup described in Section MI Al The solid 
and dashed lines correspond to the matrix elements of the 
CBF and G-matrix effective interactions, respectively, 
while the dot-dash line has been obtained replacing V c e 
with the bare v' 6 potential. Renormalization of the NN in- 
teraction, carried out either solving the Bethc-Goldstonc 
equation or modifying the basis states, appears to have 
a strong impact on the matrix elements, which become 
more attractive with respect to the matrix elements of the 
bare interaction. On the other hand, CBF and G-matrix 
approaches yield rather similar results in the considered 
range of densities and momenta. 

The matrix elements M( 3 P 2 ) = M\l are shown in Fig. 
[2J It appears that in this channel the results obtained 
using the CBF are closer to those corresponding to the 
bare interaction, while being appreciably different from 
the G-matrix results. Note, however, that the matrix 
elements corresponding to the z Pi channel are over one 
order of magnitude smaller than those corresponding to 
the 1 So channel. 

It has to be pointed out that the density dependence of 
the matrix elements is rather mild. In the CBF approach 
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FIG. 1: (Color online) Top panel: M( 1 5o) matrix element 
of the CBF (solid lines) and G-matrix (dashed lines) effective 
interactions (Moo of Eq. (|22|) 'l for neutron matter at nuclear 
matter equilibrium density, as a function of relative momen- 
tum. For comparison, the dot-dash line shows the result cor- 
responding to the bare v' e potential. Bottom panel: same as 
the top panel, but for density p = 2po. 



the density dependence arises from the two-body corre- 
lation functions, while in the G-matrix it comes through 
the presence of the Pauli operator and the single particle 
potentials appearing in the denominator of Eq.©. In 
addition, one should also take into account the density 
dependence associated with the starting energy E, as the 
matrix elements reported in the Figs. [T] and [5] have been 
computed at twice the Fermi energy of the corresponding 
density. 

The convergence of the partial wave expansion is illus- 
trated in Fig, [3J showing the ratio 
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as a function of CM energy (see Eg . (|2"0"]) ) . In the above 
equation <7 tot is the total in-medium neutron-neutron 
cross section, while a^ ot denotes the contribution of the 
£-th partial wave. All cross sections have been evaluated 
in the kinematical set up relevant to the calculation of 
the transport coefficients, discussed in Section llll Al The 
definition obviously implies that, as L — > oo, Rl — ► 1. 

The results of Fig. [3] show that the total cross sections 
obtained including only the partial waves with I = 
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FIG. 2: (Color online) Same as in Fig. Q] but for the M( 3 P 2 ) 
matrix elements. 
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FIG. 3: (Color online) Energy dependence of the ratio Rl, 
defined by Eq. (|23[) . The dashed, dot-dash and solid lines have 
been obtained including the contributions of states of angluar 
momentum £ up to 0, 1 and 2. The top and bottom pan- 
els correspond to CBF and G-matrix effective interactions, 
respectively. 



IV. TRANSPORT COEFFICIENTS 



and 1, corresponding to the spin-singlet and spin-triplet 
states of lowest angular momentum, is within less than 
5% of the fully converged result. Comparison between 
the top and bottom panels, corresponding to CBF and 
G-matrix effective interactions, also shows that the two 
approaches lead to a similar qualitative behavior, al- 
though the G-matrix Rq exhibit a somewhat steeper en- 
ergy dependence. 



A. Abrikosov-Khalatnikov formalism 

The theoretical description of transport properties of 
normal Fermi liquids is based on Landau theory Q|. 
Working within this framework and including the lead- 
ing term in the low-temperature expansion, Abrikosov 
and Khalatnikov 0] obtained approximate expressions 
for the shear viscosity and the thermal conductivity. Let 
us consider viscosity, as an example. The AK result reads 



In Fig. |4]we compare the total neutron-neutron cross 
section at nuclear matter equilibrium density, computed 
in the kinematics of Section UlI A[ as a function of Ecm- 
The difference between the cross sections obtained us- 
ing the CBF (solid line) and G-matrix (dashed line) ap- 
proaches does not exceed ~ 20 % at Eqm > 10 MeV. On 
the other hand, the screening effect due to the presence of 
the nuclear medium, illustrated by the difference between 
the solid and dashed lines and the dot-dashed one, cor- 
responding to the frcc-spacc cross section obtained from 
the i-matrix associated with the v' 6 potential, turns out 
to be large. At E C m > 100 MeV, where the CBF and 
G-matrix results are very close to one another, the in 
medium cross section turns out to be quenched by a fac- 
tor ~ 3. 



VAK = -pm v F r ^ 2(i _ ^ , 



(24) 



where Vp = kp/m* is the Fermi velocity and m* and 
t denote the quasiparticle effective mass and lifetime, 
respectively. The latter can be written in terms of 
the angle-averaged scattering probability, (W) , with (see 
Eq.©) 
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FIG. 4: (Color online) Total in-medium neutron-neutron 
cross section in neutron matter, computed at nuclear matter 
equilibrium density as a function of energy in the center of 
mass frame. The solid and dashed lines have been obtained 
using the CBF and G-matrix effective interactions, respec- 
tively, in the kinematical setup discussed in Section III A. 
For comparison, the dot-dash line shows the free-space cross 
section, obtained from the t-matrix associated with the v'q 
potential. 



where T is the temperature and 

m= [*i W 



2ir cos (6/2) ' (2?) 

Note that, as pointed out in Section UlIAl the scattering 
process involves quasiparticles on the Fermi surface. As 
a consequence, for any given density p, W depends only 
on the angular variables 6 and <f>. Finally, the quantity 
A,, appearing in Eq. (|24]) is defined as 



A r) — 



(W[l -3 sin 4 (6/2) t 
W) 



(28) 



The exact solution of the equation derived in Ref . [|| , 
obtained by Brooker and Sykes [2l[, reads 
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the size of the correction with respect to the result of 
Eq.pgj) being 0.750 < (r]/r] AK ) < 0.925. 

The expression of the thermal conductivity, k, is ob- 
tained following the same procedure, the only difference 
being that in this case the leading term in the low energy 
expansion is linear, rather than quadratic, in the inverse 
temperature T . The resulting expression is 

3-A K 

K = K AK 
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where 



and 
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(32) 



In this case the correction to the AK result turns out to 
be larger. From Eq. ([32|) it follows that — 1 < A K < 3, im- 
plying in turn (sec Eq. (00])) 0.417 < (k/kak) < 0.561. 



B. Results 

Figures[5]and[6]show the T-indcpendent quantities r]T 2 
and kT, respectively, as a function of density. The cal- 
culations have been carried out using the formalism de- 
scribed in the previous Section and the scattering prob- 
abilities W(9, </>) obtained from both the G-matrix and 
CBF effective interactions, which have been computed 
at zero temperature. For comparison, in Fig. [5] we also 
display, by the dot-dash line, the results obtained from 
the free-space scattering probability, computed using the 
t-matrix associated with the bare v' e potential. 

As the ratios m* /m resulting from the two approaches, 
CBF and G-matrix, turn out to be rather close to one 
another (the difference never exceeds few percent in the 
density range shown in Figs. [5] and [6|) all calculations 
have been carried out using the CBF effective masses. 

Figure [5] clearly indicates that medium modifications 
of the NN scattering cross sections play a critical role, 
leading to a dramatic enhancement of the viscosity. A 
similar effect is observed in the case of thermal conduc- 
tivity [22|. 
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FIG. 5: (Color online) Density dependence of the T- 
independent quantity rjT in pure neutron matter. The solid 
and dashed lines correspond to the CBF and G-matrix effec- 
tive interactions, respectively, while the dot-dash line shows 
the results obtained from the f-matrix associated with the v' 6 
potential. 
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Comparison between Fig. [5] and Fig. [6] shows that, 
while the density dependence of the thermal conductiv- 
ity resulting from the two approaches looks remarkably 
similar, in the case of viscosity sizable discrepancies oc- 
cur at densities larger than nuclear matter equilibrium 
density. 
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FIG. 6: (Color online) Density dependence of the T- 
independent quantity nT in pure neutron matter. The solid 
and dashed lines have been obtained using the CBF and G- 
matrix effective interactions, respectively. 

This feature can be easily explained considering 
the different angular dependence of the integrands in 
Eqs. (j2"5)l and (f3"2"|) , which determine the functions X v and 
X K , respectively. In the case of X K , W(0,<fr) is weighted 
with a function that only depends on 0. As a conse- 
quence, the result only depends on the ^-integrated scat- 
tering probability, which is trivially related to the to- 
tal cross section (see Eq.lpTjl) at energy Eqm given by 
Eg . (|2T))) . Hence, the similar behavior exhibited by the 
two curves of Fig. [6] merely reflects the fact that the to- 
tal cross sections obtained from the G-matrix and CBF 
approaches turn out to be close to one another (see Fig. 
[2J . On the other hand, in the right hand side of Eq. ([25]) 
the scattering probability is weighted with a factor that 
emphasizes the differences in the <j) dependence of the 
CBF and G-matrix W(6, (f>). 

We have verified that the large discrepancy between 
the solid and dashed curves of Fig. \5\ at high density 
is in fact ascribable to X v . The A^-independent quanti- 
ties t]ak{^ — X v )/4 obtained from the CBF and G-matrix 
approaches turn out to be within less than 5 % of one 
another at p = 0.32 fm -3 . On the other hand, removal 
of the X K dependence in kak does not produce any sig- 
nificant effects. 



V. CONCLUSIONS 

Many body theory provides a fully consistent frame- 
work, suited to construct effective interactions starting 



from highly realistic models of nuclear dynamics. In 
this work, we have employed the effective interactions re- 
sulting from the G-matrix and CBF approaches to com- 
pute the in-medium NN cross sections, which are needed 
to obtain the transport coefficient within the Landau- 
Abrikosov-Khalatnikov formalism. 

The calculations have been carried out using the trun- 
cated v' 6 form of the NN potential of Ref . [ll[ ■ The effects 
of three- and many-body forces, though being known to 
be sizable at large density, have not been taken into ac- 
count. The results of Ref. show that many-body 
forces give rise to a change of the shear viscosity of less 
than 10 % at p < 0.32 fm~ 3 . Hence, their inclusion is 
not likely to significantly affect the main conclusions of 
our work. 

The approach based on effective interactions allows 
one to consistently take into account screening effects 
arising from short range NN correlations, that lead to a 
large decrease of the scattering cross section. As a con- 
sequence, the shear viscosity and thermal conductivity 
obtained from the effective interactions turn out to be 
much larger than the corresponding quantities computed 
using the bare NN potential. 

Our work, showing that the results of the G-matrix 
and CBF schemes are in reasonable agreement with one 
another, suggests that as long as the effective interaction 
is based on a realistic NN potential, strongly constrained 
by the large data set of NN scattering phase shifts, the 
model dependence associated with the many body ap- 
proach employed is not critical. 

On the other hand, it has to be pointed out that 
Skyrme effective interactions (for a recent and compre- 
hensive discussion the application of the Skyrme ap- 
proach to nuclear matter and neutron stars see Ref. p3| ) . 
mainly constructed by fitting bulk properties of nuclear 
matter, predict in-medium NN cross sections whose be- 
havior is significantly different from the one predicted by 
the G-matrix and CBF effective interactions. As a re- 
sult, the values of the viscosity and thermal conductivity 
coefficients computed using Skyrme effective interactions 
turn out to be much lower than those shown in Figs. [5] 
and [SI For example, using the SLya effective interaction, 
adjusted to reproduce the the microscopically derived 
EOS of neutron and nuclear matter (24|, one finds at nu- 
clear matter equilibrium density r/T 2 ~6x 10 13 g cm -1 
s _1 MeV 2 and kT ~ 4 x 10 30 erg cm -1 s -1 , to be com- 
pared to ~ 1.4xl0 15 gem" 1 s _1 MeV 2 and kT ~ 4xl0 31 
erg cm -1 s _1 obtained from the G-matrix and CBF for- 
malisms. While the Skyrme approach has proved to be 
very useful in many contexts, these results suggest that 
the determination of the transport properties of nuclear 
matter requires effective interactions providing a quanti- 
tative account of the observed NN scattering data in the 
limit of vanishing density. 
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